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ABSTRACT 

We report on atomic gas (H I) and molecular gas (as traced by CO (2-1)) redshifted absorption 
features toward the nuclear regions of the closest powerful radio galaxy, Centaurus A (NGC 5128). 
Our H I observations using the Very Long Baseline Array allow us to discern with unprecedented 
sub-parsec resolution H I absorption profiles toward different positions along the 21 cm continuum jet 
emission in the inner O'.'S (or 5.4 pc). In addition, our C0(2-l) data obtained with the Submillimeter 
Array probe the bulk of the absorbing molecular gas with little contamination by emission, not possible 
with previous CO single-dish observations. We shed light with these data on the physical properties 
of the gas in the line of sight, emphasizing the still open debate about the nature of the gas that 
produces the broad absorption line (^ 55 km s~^). First, the broad H I line is more prominent toward 
the central and brightest 21 cm continuum component than toward a region along the jet at a distance 
^ 20 mas (or 0.4 pc) further from it. This suggests that the broad absorption line arises from gas 
located close to the nucleus, rather than from diffuse and more distant gas. Second, the different 
velocity components detected in the CO (2-1) absorption spectrum match well other molecular lines, 
such as those of HCO+(1-0), except the broad absorption line that is detected in HCO"''(1-0) (and 
most likely related to that of the H I) . Dissociation of molecular hydrogen due to the AGN seems to be 
efficient at distances r < 10 pc, which might contribute to the depth of the broad H I and molecular 
lines. 

Subject headings: galaxies: elliptical and lenticulars — galaxies: individual (NGC 5128) — galaxies: 
structure — galaxies: ISM 



o 

O 

o 



% 



1. INTRODUCTION 

The analysis of absorption lines arising from atomic 
and molecular gas in galaxies is a powerful technique to 
study the physical and chemical properties of the inter- 
stellar medium (ISM) . Only a handful of nearby galaxies 
are known to show atomic (H I) and molecular absorp- 
tion features. H I and molecular absorption line studies 
are not only relevant to deduce the properties of the ISM 
along specific lines of sight in nearby galaxies (e.g. Gupta 
et al. 2007), but are also essential for an interpretation of 
similar data toward higher redshift galaxies (e.g. Muller 
et al. 2006; Combes 2008; Henkel et al. 2009). Studying 
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both H I and molecular gas absorption lines at the same 
time can reveal the interface between both components 
of the ISM at the circumnuclear regions where dissocia- 
tion of molecular gas is expected as a result of energetic 
processes triggered by Active Galactic Nuclei (AGNs). 

At a distance of D ~ 3.8 Mpc ^ (Harris et al. 2009), 
Centaurus A (Gen A, NGC 5128) is the nearest giant 
elliptical exhibiting nuclear activity, with powerful radio 
jets extending almost 10° on the sky and with prominent 
atomic and molecular emission as well as absorption fea- 
tures detected toward its nuclear regions (Israel 1998). It 
is thus a unique source to study the molecular to atomic 
interface. At this distance, high spatial resolution can be 
achieved: 1" corresponds to 18 pc projected linear size. 
For a comprehensive general review of Cen A, see Israel 
(1998). 

One of the most prominent features in Cen A is the 
dust lane along its minor axis, which is composed of a 
large amount of gas and dust within a warped disk-like 
structure seen nearly edge-on. An H I mass of (10 ± 
3) X 10^ M0 extends out to a radius of at least r — 
7' (or about 7 kpc), while it is observed to be mostly 
absent within the inner arcmin (or 1 kpc projected linear 
size) (van Gorkom et al. 1990; Schiminovich et al. 1994; 
Israel 1998). A comparable amount of molecular gas is 
contained within the inner r ~ 2' (2 kpc), Mh2 — 4.2 
xlO^ Mq (Eckart et al. 1990b). About 30% - 40% of the 
total molecular gas content is located within the inner 1' 
(1 kpc), Mh2 = (8.1 ± 0.4) xlO^ Mq (Espadaet al. 2009, 
hereafter Paper I) . A circumnuclear disk of molecular gas 

1 This distance is used throughout this paper. 
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is resolved in the inner kiloparsec, with a size of 24" x 12" 
(430 X 215 pc) in diameter, and with different kinematics 
with respect to the molecular gas at larger radii (Paper 

I). 

Along the NE - SW direction (perpendicular to the 
circumnuclear molecular gas disk) , there exists a spectac- 
ular radio-continuum jet extending from the nucleus (pc 
scale) to the outer radio lobes (hundreds of kpc scale). 
Subparsec-scale VLBI (Very Large Baseline Interferom- 
etry, including the VSOP - VLBI space observatory pro- 
gramme) high resolution radio-continuum observations 
(e.g. at 2.3, 4.8, 8.4, 22.2 and 43 GHz) show that the 
structure of Cen A is rather complex, consisting of a 
bright nuclear jet and a fainter counterjet, as well as a 
bright central component (e.g. Tingay et al. 1998; Tin- 
gay & Murphy 2001; Horiuchi et al. 2006). Note that this 
component is not the nucleus itself since the nucleus is 
not bright at frequencies below 10 - 20 GHz due to syn- 
chrotron self-absorption and free-free absorption from a 
disk or torus of ionized gas (Jones et al. 1996; Tingay 
& Murphy 2001). The compact core becomes visible at 
mm/submm wavelengths (Kellermann et al. 1997) where 
its spectral index is quite flat (see Fig. 6 in Meisenheimer 
et al. 2007). 

H I absorption lines in this source were found for the 
first time four decades ago (e.g. Roberts 1970). Since 
then, VLA aperture synthesis observations with resolu- 
tions of the order of ~ 10" (van der Hulst et al. 1983; 
Sarma et al. 2002) have shown H I absorption lines in 
different regions along the jet. They are distributed to- 
ward both the core and the inner jet (extended along 30") 
and are composed of a set of lines with widths of typi- 
cally 10 km s^^. The most prominent H I absorption 
component is located close to the systemic velocity of 
the galaxy (Vlsr = 541.6 ± 2 km s-\ Israel 1998)^. To- 
ward the compact bright continuum component, two red- 
shifted lines are seen at V^ = 578 km s~^ and 598 km s^^, 
as well as a broad absorption component (extending up 
toV = 620 km s^^ with a width of AF ~ 55 km s^^). 
Sarma et al. (2002) showed that the latter is only present 
toward the central component, but not along the other 
regions of the jet. This sets an upper limit to the size of 
the absorbing material of about 20", or 360 pc as seen in 
projection. 

Molecular absorption lines from cm to submm wave- 
lengths have also been detected for several species. At 
cm wavelengths narrow molecular lines (a few km s^^ 
width) close to the systemic velocity have been observed 
in several molecules such as OH, H2CO, C3H2 and NH3 
(Gardner & Whiteoak 1976a,b; Seaquist & Bell 1990; van 
Langevelde et al. 1995, 2005). Single-dish observations 
have been performed at mm/submm wavelengths for the 
lowest transitions of CO, HON, HNC, HCO+, Hi3C0+, 
N2H+, C2H, CN and CS (Eckart et al. 1990a; Israel et al. 
1990; Wiklind & Combes 1997; Eckart et al. 1999). In 
addition to narrow lines close to the systemic velocity, 
both narrow and broad absorption features are seen at 
redshifted velocities in HCO+(1-0), OH 18cm, and to 
a minor extent in HCN(l-O), HNC(l-O) and CS(2-1), 
covering a velocity range from ^ 570 to 620 km s~^ (e.g. 
Seaquist & Bell 1990; Eckart et al. 1990b; Israel et al. 

^ All velocities in this paper are radio LSR. Conversion from 
LSR to heliocentric velocity is VhcI = ^LSR + 2.4 km s^^. 



1991; Wiklind & Combes 1997). 

Although there is a consensus that narrow absorption 
lines with velocities close to the systemic velocity repre- 
sent gas far from the nucleus, a dichotomy exists regard- 
ing the nature of the broad redshifted absorption line. It 
has been suggested that it might be: i) warm and dense 
gas close to the nucleus and probably falling into it (e.g. 
van der Hulst et al. 1983; Israel et al. 1991), and ii) cold 
diffuse gas above and below a warped disk (Eckart et al. 
1999). In this paper we aim to shed light on the model 
that better represents the different gaseous components 
(both atomic and molecular) producing such absorption 
lines, and their physical properties as a result of their 
location with respect to the nucleus. In particular, we 
emphasize the origin of the broad line component. 

We have carried out H I observations using the Very 
Large Baseline Array (VLBA) and the phased Very Large 
Array (VLA)^ which allow us to derive the properties of 
the main H I absorption features seen toward the con- 
tinuum source at 21 cm with unprecedented sub-parsec 
resolution (versus the ^ 180 pc resolution in past VLA 
experiments). This is complemented by C0(2-l) obser- 
vations (and isotopologues) at the Submillimeter Array 
(SMA'', Ho et al. 2004) (Paper I) from which we can 
extract a C0(2-l) absorption profile with minimal con- 
tamination by emission. We introduce our VLBA H I 
and SMA C0(2-l) observations and the data reduction 
in § 2, where we present the continuum maps and ab- 
sorption spectra. In § 3 we focus on the identification of 
the individual H I absorption components and the physi- 
cal properties of their corresponding regions. We present 
our results on the C0(2-l) and ^^C0(2-l) absorption 
lines in § 4 and compare the H I absorption lines with 
those of molecular lines in § 5. We discuss the results in 
§ 6 and summarize our findings in § 7. 

2. OBSERVATIONS AND DATA REDUCTION 

The array configuration, phase center, primary and 
synthesized beam, achieved noise and spectral resolution 
of our VLBA H I and SMA C0(2-l) observations are 
shown in Table 1. In the following we explain in more 
detail our observational strategies and data reduction. 

2.1. VLBA HI Observations 

The VLBA H I observations were carried out on 26 
November 1994. Due to its low declination (Dec. ~ 
—43°), Cen A was observable for 2 hours with 8 of the 10 
VLBA antennas and the phased VLA. The observations 
were made using 2 IPs of 2 MHz each. The IPs were 
centered on the redshifted H I line (t'rost = 1420.405752 
MHz, at a velocity Vlsr — 541.6 km s^^). Global fringe- 
fitting was performed using AIPS^ task PRING with a 
solution interval of 6 minutes. Amplitude calibration was 
derived using the antenna gain and system temperature 
recorded at each antenna. Bandpass calibration was done 

■^ The VLBA and VLA are operated by the National Radio As- 
tronomy Observatory, a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 

■* The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

^ http://www.aips.nrao.edu/cook.html. 
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using the calibrator 3C273 (J1229+0203). The contin- 
uum was obtained from line-free channels. Continuum 
subtraction was done in Difmap (Shepherd et al. 1995) 
using a model of CLEAN components made from several 
line-free channels. Subsequent editing and imaging of all 
data was also done using Difmap, using natural weight- 
ing. The resulting Half Power Beam Width (HPBW) is 
36.3 X 7.2 mas (0.65 x 0.13 pc). The velocity resolution 
in these observations is 3.3 km s~^, and the total velocity 
coverage of a single IF is ^ 430 km s~^. 

2.2. SMA C0(2-l) Observations 

Details of the C0(2-l) observations using the SMA in 
its compact configuration are reported in Paper I. C0(2- 
1), i^C0(2-l) and 0^^0(2-1) lines were simuUaneously 
observed. We used two independent sources, 3C273 and 
Callisto, to check the consistency of our bandpass cali- 
bration. Both solutions were nearly identical in most of 
the spectral windows, and in particular where the C0(2- 
1) lines were located. The continuum toward CcnA at 
1.3 mm was found to consist of an unresolved source with 
a flux of 5*1. 3mm = 5.9 ± 1.0 Jy beam~^ at the AGN po- 
sition which was used to calibrate the data of Cen A it- 
self (Paper I). No continuum subtraction was performed. 
The velocity resolution is 1 km s~^. 

In order to minimize undesired emission, we obtain the 
CO (2-1) absorption spectrum by directly fitting the am- 
plitude of the interferometric visibilities using a point 
source model fixed at the position of the unresolved con- 
tinuum source. This was done with the task UV_FIT in 
the GILDAS^ reduction package. The C0(2-l) spectrum 
extracted using visibilities from all baselines still shows 
some emission at velocities V = 420 - 510 km s^^ and V 
= 620 - 680 km s~^. This corresponds to extended emis- 
sion from the circumnuclear regions (Paper I). A new fit 
to the visibility data with projected baseline lengths > 
50 m yielded a flatter spectral baseline. This is due to 
a combination of the spatial flltering capabilities of the 
interferometer, beam smearing and probably a lack of 
C0(2-l) emission toward the central ^ 5" [^ 90 pc). 
The typical extent that we take into account by remov- 
ing the shorter baselines can be roughly estimated as 
< \/D ~ 4" [^ 72 pc), where A is the wavelength of the 
observations and D the average baseline length. 

We show in Figure 1 the resulting C0(2-l) spectra 
with and without visibilities with baseline lengths larger 
than 50 m. For the 13C0(2-1) and 0^^0(2-1) transi- 
tions, however, we retained all visibilities since there is 
no signiflcant emission down to our noise level. Their 
resulting noise per channel, cr = 0.1 Jy, is thus smaller 
than in the long-baseline only C0(2-l) spectrum {a — 
0.3 Jy). 

3. HI ABSORPTION FEATURES 

3.1. 21 cm Continuum: Nuclear Jet 

The 21 cm continuum map in the inner 0'.'3 (5.4 pc) 
is presented in Figure 2 (center). This high resolution 
map resolves the bright 21 cm central component {^ 
10", or 180 pc) found by van der Hulst et al. (1983) 
and Sarma et al. (2002) (labeled as "Nucleus"). We flnd 
that this central component further comprises a more 

6 http://www.iram.fr/IRAMFR/GILDAS. 



compact bright component and a nuclear jet extending 
up to 0'.'23 to the NE, with a simflar P.A. ~ 51° as the 
inner 21 cm jet extending up to 55" (van der Hulst et al. 
1983; Sarma et al. 2002). Note that van Langevelde et al. 
(2005) present a VLBI 18 cm continuum map with a sub- 
arcsecond resolution where the central bright component 
is seen. However, the rest of the nuclear jet is not as 
clearly detected as in our map. 

While the peak flux density is 603 mJy beani^^ and 
the noise level of our maps is 2 mJy beam~^, note that 
the uncertainty of the absolute fiux density measurement 
is larger. Uncertainties in flux calibration are estimated 
to be about 5% (5*2 icm = 600 ± 30 mJy beam~^). In 
addition, on-source errors due to deconvolution are ex- 
pected. For instance, Tingay & Murphy (2001) estimated 
that the corresponding flux density error in their 2.2 GHz 
VLBA data of Cen A due to this effect is about 30% for 
the brightest components. This source of error is ex- 
pected to be larger for our lower observing frequency, 
since the uv coverage is less uniform. 

As mentioned in § 1, the core of the radio source has 
a very strongly inverted spectrum since the continuum 
emission is seen through a circumnuclear ionized gas disk 
which is opaque at wavelengths longer than 13 cm (or fre- 
quencies lower than 2.3 GHz) due to free-free absorption 
and to self-absorbed synchrotron emission (Jones et al. 
1996; Tingay & Murphy 2001). The compact bright com- 
ponent in Figure 2 is thus part of the approaching jet 
and is offset from the actual AGN. Tingay & Murphy 
(2001) showed that there is an absolute offset of 12.5 mas 
to the E and 10.5 mas to the N between the 8.4 GHz 
and 2.2 GHz central components. This corresponds to 
a separation of 0.3 pc along the jet axis. We show in 
Fig. 2 the estimated location of the AGN (8.4 GHz cen- 
tral component), assuming that the bright component of 
our 1.4 GHz map coincides with that of the 2.2 GHz map, 
and also that position variability between the different 
observing dates is not significant. We do not detect any 
continuum emission from the counter-jet toward the SW 
to the 2 mJy beam"-'^ noise level of our VLBA maps. 

3.2. Identification of H I Absorption Features 

Figure 3 shows the H I spectrum integrated over the 
entire continuum source. We find 5 main H I absorp- 
tion lines: four relatively narrow lines and an underlying 
broad line. These absorption lines, all red-shifted with 
respect to the systemic velocity, were detected using the 
VLA by van der Hulst et al. (1983) and Sarma et al. 
(2002), although toward an unresolved component of size 
~ 8" - 10". We adopt the nomenclature used in Sarma 
et al. (2002) for the most prominent lines (in ronian nu- 
merals, IV, SV, I, II and HI). 

We derive from Gaussian fits the mean velocities (V^), 
peak flux densities {S) and Full Widths at Half Maxi- 
mum (FWHM) of the different lines in our H I proflle, 
as shown in Table 2. We first fit the broad line using 
the redshifted wing &tV > 600 km s~^, and then fit the 
remaining lines individually. An exception is line I/H, 
which is a blend of lines I and II in the H I spectra of 
Sarma et al. (2002). Being almost unresolved with our 
3.3 km s~^ resolution, we decided to compute them to- 
gether. The H I absorption features in our spectra are 
located sXV = 541, 552, 575 and 595 (± 1) km s^^ (IV, 
SV, I/H, HI, as outlined in Figure 3). The most promi- 
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nent lines show FHWM values of about 5-10 km s ^. 
A fit to the broad absorption line yields a FWHM = 53 
± 10 km s~i centered at 578 ± 4 km s"^ The final fit 
and the resulting residual from the observed H I profile 
are presented in Figure 3. 

The great advantage of our VLBA data is that we 
can obtain the H I absorption spectrum toward differ- 
ent positions over the continuum emission. Figure 2 
shows the four most prominent continuum sources where 
we obtained H I absorption profiles. Their relative off- 
sets with respect to the phase center are: {Aa,A6) — 
(+73,+69) mas, (+63,+58) mas and (+20,+10) mas 
(namely, positions a, b and c), corresponding to sep- 
arations of 1.8 pc, 1.5 pc and 0.4 pc respectively; as 
well as toward the brightest component at (-1-1,-1) mas 
(namely, position d). The continuum flux density peaks 
of the different components are 52icm = H, 23, 145 and 
603 mjy beam~^, respectively. 

While we detect H I absorption features only close to 
the systemic velocity (components IV and SV) in posi- 
tions a and b, in position c and d the profiles show two 
additional redshifted H I absorption components I/II and 
III at F = 574 - 576 km s'^ and 595 - 596 km s~^ (Fig- 
ure 4). An additional shallow redshifted wing in position 
d is clearly seen at velocities V > 600 km s~^, which cor- 
responds to the broad absorption line. The fitted Gaus- 
sian parameters are presented in Table 3, and the fits 
are shown in Figure 4. To fit these lines, we set as initial 
values the velocities and FWHM found for the H I ab- 
sorption lines integrated over the entire continuum. Since 
the flux densities refer to their corresponding continuum 
emission, we calculated next the optical depths and H I 
column densities in order to compare between different 
positions. 

3.3. H I Optical Depths and Column Densities 

We calculated the H I optical depths, thi, as Iohs[Jy] — 
IcontS~^™, where Jobs is the observed fiux in absorption 
and /cont is the measured continuum flux. Absolute flux 
uncertainties do not have any effect on the optical depth 
since it is derived from the /obs/^cont ratio. In columns 5 
and 6 of Table 2 and 3 we present the H I optical depth 
peak (peak thi) and integrated optical depths {J TindV) 
integrated over the entire continuum source and for each 
independent position, respectively. The upper limits cor- 
responding to non-detected lines are calculated as 3 cr 
{AV Sv)^/"^, where AV is the FWHM in Table 2 and dv 
is the channel spacing of 3.3 km s^^. 

In Figure 5 we plot thi toward the four different posi- 
tions a, 6, c and d along the 21 cm nuclear jet. The main 
absorption line is component SV, with peak thi > 0.75 
seen at all positions. No trend is seen as a function of 
distance from the nucleus. On the other hand, the red- 
shifted velocity component I/II seems to be preferentially 
deeper if closer to the brightest component (position d). 
Component I/II is characterized by peak thi = 0.29 ± 
0.03 in position c, while peak thi — 0.41 ± 0.02 in posi- 
tion d. The 3(7 upper limit obtained in position b is peak 
Thi < 0.3. However, we would have expected at least a 
tentative detection for component I/II in position b if of 
comparable depth as in position c or d. 

It is surprising that we detect the broad line compo- 
nent against the position d but it is weaker toward the 
adjacent continuum region, position c. Figure 6 shows 



a close-up of the thi in both locations. A fit to the 
broad component (assuming that a Gaussian function 
applies) gives integrated optical depths of / TnidT^ — 3.6 
± 1.3 km s^^ for position d while it is J TuidV — 0.5 ± 
1.1 km s~^ in position c (Table 3). That the broad ab- 
sorption line in position c is weaker is reinforced by the 
smaller optical depths in the spectral regions between the 
hues SV, I/II and HI atV ^ 560 km s^^ and 590 km s"! 
with respect to position d (Figure 6). This implies that 
the H I opacity of the broad line component varies within 
a projected linear scale of less than 0.4 pc. 

Finally, we calculate the H I column densities (Column 
7 in Tables 2 and 3), using the following relation: 



7V(HI)/T, [cm"2 j^- 



1.823 X 10 



18 



I midV [kms-i] 



where Tg is the spin temperature. The H I column densi- 
ties in the different lines cover the range 7V(H I)/Ts ~ (1 
- 12) X 10^* cm~^ K~^. The broad line toward position 
d has a corresponding A^(H I)/Ts that is comparable to 
the deepest narrow lines. 

4. CO (2-1) ABSORPTION FEATURES 

In our 1.3 mm observations, with a restoring beam of 
6'.'0 X 2'.'4 (or 110 x 40 pc), we found an unresolved con- 
tinuum source at the galactic center (Paper I). The inner 
jet observed at cm wavelengths has a steep radio spec- 
trum so the contribution at mm wavelengths is negligi- 
ble (Israel et al. 2008). Therefore our 1.3 mm continuum 
emission is dominated by the nuclear continuum emis- 
sion. The angular size of the central continuum source 
is 0.5 ± 0.1 mas at 43 GHz, or 0.009 ± 0.002 pc (Kcller- 
mann et al. 1997), and is considerably smaller than our 
synthesized beam. The absorption spectra of C0(2-l), 
i3CO(2-l) and 0^^0(2-1) (the latter two taking into ac- 
count all visibilities) are shown in Figure 7 6), c) and d), 
respectively. Absorption features have been detected in 
C0(2-l) and 13C0(2-1), but not in Ci*0(2-1) at a 3 ct 
level (= 0.3 Jy). 

4.1. C0(2-l) Line Characterization 

The C0(2-l) absorption spectrum spans a range of 
540 < T^ < 620 km s~^ and is composed of 6 narrow 
lines. The narrow absorption features are labelled as 
outlined in Figure 7 a): lines 1, 2 and 3 correspond to 
the absorption lines close to the systemic velocity, at V 
= 539, 544 and 550 km s~^, and 4^ 5 and 6 to absorption 
lines at higher velocities, located &t V = 575, 606 and 
613 km s^^ (Table 4). Component 4 is clearly detected 
in the profile with all visibilities, but it turns out to be 
weak with only the longest baselines due to the increase 
of noise level. C0(2-l) absorption line 2 is located in the 
most likely systemic velocity of the galaxy. Absorption 
lines i, 2 and 3 are much deeper than the higher velocity 
components, with component 3 the most prominent. 

Three ^^C0(2-l) absorption lines are also detected, 
corresponding to features 1, 2, 3 in the C0(2-l) profile. 
An additional line corresponding to C0(2-l) line 5 was 
tentatively detected. Features 4 a-nd 6 are not detected 
at a 3 (7 level. 

The linewidths are comparable in most cases to our 
spectral resolution. Because of the limited velocity reso- 
lution of our SMA observations (1 km s^^) and the nar- 
rowness of the absorption lines, we do not attempt to fit 
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Gaussian functions to the several C0(2-l) Hne compo- 
nents. Line widths in the C0(2-l) profile of components 

1 and 3 are about Ay = 3-5 km s~^, while the re- 
maining lines have even smaller line widths AV ~ 1 - 

2 km s^^. Lines with widths less than 1-2 km s~^ are 
not (or insufficiently) resolved. Note that in those cases 
the optical depths are underestimated. 

4.2. C0(2-l) Optical Depths and Column Densities 

Since the radio core is small {^ 0.009 pc, KcUcrmann 
et al. 1997), the covering factor is likely close to unity. In 
fact, the components close to the systemic velocity (i, 2 
and 3) have depths > 0.8, and assuming that / > depth 
of absorption for a normalized spectrum (continuum flux 
is unity), the covering factor is / > 80%. For simplicity 
we adopt a filling factor / of unity for all the components. 

We have calculated the C0(2-l) optical depths 
'^co(2-i), in the same manner, Iohs[Jy\ = ^conte~'^<=o(^-i) 
(Table 4), as those of H I. The optical depths for any 
undetected "C0(2-l) and 0^^0(2-1) lines have been 
calculated using a 3 cr level upper limit. C0(2-l) com- 
ponent 3 is saturated since the depth of the absorption 
line nearly equals the continuum level, and tco(2-i) is 
highly uncertain due to non-linearity of the logarithmic 
function. High values of peak Tco(2-i) show that com- 
ponents 1, 2 and 3 are optically thick. C0(2-l) compo- 
nents 4, 5 and 6, as well as all of the "'^^00(2-1) lines, are 
characterized by peak tco(2-i) ^ 0.8. Thus, they can be 
considered as optically thin if filling factor is unity. 

Our C0(2-l) optical depths arc reasonably well in 
agreement with previous single-dish measurements. The 
peak rco(2-i) for components i, 2 and 3 are tco(2-i) 
= 2.1 ± 0.3, 1.6 ± 0.2 and > 3.6 ± 1.6, respectively, 
which are compatible with those obtained by Israel et al. 
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(1991), rco(2-: 
(1991) also gives an optical depth for component 4 of 
Tco(2-i) = 0.4 ± 0.2, while our value is 0.10 ± 0.05. The 
remaining lines were not detected by Israel et al. (1991). 
This comparison shows the need of interferometric obser- 
vations in order to minimize confusion by emission and 
obtain higher accuracy. 

As for the ^'^00(2-1) lines, the peak ti3co(2-i) values 
are: 0.11 ± 0.02, 0.19 ± 0.02, 0.58 ± 0.03 and 0.05 ± 
0.01 for components 1, 2, 3 and 5. Eckart et al. (1990a) 
gives a peak optical depth of 0.4 for the line close to F = 
550 km s~^ and non-detections for the rest of the lines. 
This value is well in agreement with our component 3. 

We derive column densities of each C0(2-l), ^■^C0(2- 
1) and C^^0(2-l) lines using the following equation: 

AT(rn^ 8^^3 g(fcre.)exp(i;j/fcTe.) [ 
"^^^^^ ^ c3A.,5« (l~exp(-WfcTe.)) j ^-o(-i)^^' 
where we assume that LTE conditions apply. Ej de- 
notes the energy of the lower level of the transition, 
Tex the excitation temperature, Q{Tcx) — ^t{'^J + 
1) exp (— _Ej/fcTox) the partition function, and A^i, the 
Einstein coefficient of each isotopologue. We use Tex — 
10 K, although it could vary depending on the properties 
and location of the molecular clouds. Indeed, a gradient 
in excitation temperature has been inferred using a nom- 
inal temperature ratio of C0(2-T) and CO(l-O) emission 
line ratios by Israel et al. (1990, 1991): Tex ~ 10 K at 
r > 1' (> 1 kpc) while it seems to be larger in the cir- 
cumnuclear gas (r < 20"), Tox — 25 K. 



The different detected lines show CO column densi- 
ties N{CO) spanning two orders of magnitude (Table 4), 
from (23 ± 3) x lO^"^ cm'^ (component ^) to > 1530 
X 10^^ cm^^ (component 3). N{CO) increases approx- 
imately an order of magnitude if a Tex = 100 K is used. 
The ^'^CO column density is in the range < 8 x 10^'* 
cm^^ (component 4) and (145 ± 4) x 10^'* cm^^ (com- 
ponent 3). 

We list in Table 4 the H2 column densities (iV(H2)) 
assuming abundance ratios [C0]/[H2] — 8.5 x 10^"'' (Fr- 
erking et al. 1982), ["COJ/pa] = 1 x 10"^ (Solomon 
et al. 1979) and [Ci^OJ/pa] = 1.7 x 10"^ (Frerkinget al. 
1982). The II2 column densities as obtained from C0(2- 
1) should be similar to those obtained with ^"^00(2-1). 
However, N{}i2) obtained with C0(2-l) are systemati- 
cally smaller than those obtained using "^00(2-1). The 
[CO]/[^'^CO] ratios for the six lines in order of increas- 
ing velocity are: 17 ± 2, 5 ± 1, > 10, > 3, 17 ± 7 and 
> 5. These abundance ratios are considerably smaller 
than the assumed [CO]/[^'^CO] ~ 85, similar to the lo- 
cal value. This deviation in the abundance ratios might 
be due to a combination of the following reasons: i) dif- 
ferent excitation temperatures for ^^CO and ^^CO, ii) 
optical depth effects. Hi) different filling factors and iv) 
chemistry, such as selective photodissociation (of ^■^CO). 
However, low abundance ratios of ^ 25 have been ob- 
served towards the Galactic center and PKS 1830-211 
(see Table 7, MuUer et al. 2006) 

4.3. /12 C7//13 C] and /^^ oy//i8 Q] Isotopic Ratios 

We can roughly estimate the [^^C]/p^C] and 
[iSO]/[i^O] isotopic ratios from the [CO]/["CO] and 
[^■^CO]/[C O] ratios, assuming that the molecular abun- 
dance ratios can be decomposed into isotopic abundance 
ratios of the atoms. If we assume that [^^CO]/[^'^CO] 
= [i2C]/["C], then we obtain that [^^C]/[^^C] ~ 5 - 
17. The [^^C]/[^'^C] ratio is a measure of the primary 
to secondary processing: while ^^C mostly arises from 
massive stars, ^^C is predominantly produced from re- 
processing of -'^^C from earlier stellar generations (e.g. 
Wilson & Rood 1994). 

From our observations we derive lower limits from 
the abundance ratio [^'^CO]/[C^^O]. It ranges from >2 
to >12. This is somewhat less restrictive than the 
[i-'^CO]/[Ci^O] > 30 obtained by Wild et al. (1997) us- 
ing the J=l-0 absorption line, although their spectral 
resolution is a factor 2 lower (^ 1.9 km s^^). Higher 
spectral resolution observations, with a channel width of 
0.23 km s~^, had been previously carried out by Israel 
et al. (1991). Israel et al. report a 2 tr detection in the 
C^^O line for component 5, which gives a value for the 
abundance ratio [i3C0]/[Ci*^0] ~ 13 ± 3. This would 
be in agreement with our lower limit, but not with that 
in Wild et al. (1997). 

We can then obtain p^0]/[^^0] from our previous es- 
timate of [^^C]/[^'^C]. ^ O and ^^O form in massive stars 
and are predicted to behave similarly, except at early 
stages of massive star formation since ^^O is a primary el- 
ement and ^^O is a secondary one (Prantzos et al. 1996). 
The ratio is therefore expected to decrease with time and 
stellar processing (e.g. MuUer et al. 2006). 

By adopting [12C]/["C] =5-17 and ["CO]/[Ci^O] 
~ 13 ± 3 (Israel et al. 1991) (component 3), then 
[i^O]/[i^O] - 70 - 220 . This ratio would be in agree- 
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mcnt with values ~ 150 — 200 in nearby starburst galax- 
ies such as NGC 253 and NGC 4945 (Wang et al. 2004; 
Harrison et al. 1999; Henkel et al. 1993). 

OveraU, the [i2C]/["C] and [^^0]/[^^0] ratios in Ccn A 
are within those found toward the Galactic center and 
starburst galaxies (Henkel et al. 1993). Thus, our 
calculated [i2C]/["C] and [^'^0]/[^^0] ratios using CO 
lines suggest within the uncertainties that the molecu- 
lar clouds producing the absorption lines close to the 
systemic velocity might have been affected by relatively 
long-term starburst conditions. 

On the other hand, note that Wiklind & Combes 
(1997) derived an isotopic ratio [^2C]/[^^C] > 70 for com- 
ponent 3 from its [HCO+]/[Hi^CO+] ratio, which yields 
[i*^0]/[i^O] > 200 - 840. This seems to contradict our re- 
sults and shows the need of a systematic chemical study 
using different isotopic molecular lines. 

5. COMPARISON OF ATOMIC AND MOLECULAR 
ABSORPTION LINES 

In this section we compare the variations of 
THCO+(i-o)/'rco(2-i) and tco(2-i)/7"hi in order to infer 
the variation of molecular to atomic abundances with 
velocity. As a zeroth-order approximation we consider 
three assumptions: i) absorption lines of different species 
(both molecular and atomic) that have similar veloci- 
ties are physically related; ii) temperature (and spin) 
differences vary similarly for different species so that 
the ratios are reasonably constant; and Hi) molecular 
abundance ratios are reasonably constant ([C0]/[H2] and 
[HC0+]/[H2]). 

5.1. Molecular gas as traced by C0(2-l) and 
HCO+(1-0) 

The different velocity components detected in the 
CO (2-1) absorption spectrum correspond well to other 
molecular lines, such as those of HCO"'"(1-0). We chose 
HCO+(1-0) for comparison because it traces molecular 
gas with different physical conditions (its critical den- 
sity is two orders of magnitude larger than that traced 
by C0(2-l)). It is also the best S/N and the high- 
est spectral resolution (~ 0.2 km s^^) absorption profile 
that can be found in the literature (Wiklind & Combes 
1997). The C0(2-l) absorption profile can be compared 
to that of HCO"'"(1-0) (Figure 7 e), which results from 
the sum of 17 Gaussian components used by Wiklind & 
Combes (1997) to parametrize the observed HCO+(1-0) 
absorption spectrum. The more prominent lines of the 
HCO+(1~0) spectrum are centered at the same veloci- 
ties and have similar widths to their CO (2-1) counter- 
parts. Note that no variability was found in different 
epochs for the HCO+(1-0) absorption spectrum (Wik- 
lind & Combes 1997), so we do not expect to find vari- 
ability in our C0(2-l) profiles either. 

On the other hand, the high velocity component of 
the HCO+(1-0) absorption profile (Wiklind & Combes 
1997; Eckart et al. 1999) (Figure 7 e) is more prominent 
than in our C0(2-l) spectrum (Figure 7 h). The broad 
component of HCO"'"(1-0) absorbs roughly 15% of the 
continuum emission (Wiklind & Combes 1997). If the 
C0(2-l) broad absorption feature were as deep as the 
HCO+(1-0) broad line then we would have detected it. 

We derive the abundance ratio assuming T^x = 10 K 
for C0(2-l) and To^ = 5 K for HCO+(1-0) (Wiklind 



& Combes 1997). For fines within 540 km s"^ V < 
570 km s^^, the total abundances of HCO+ and CO 



are 7V(HC0^ 



4.8 X 10" cm-^ (Table 8, Wiklind 



& Combes 1997) and iV(CO) ~ 2.3 x 10^^ cm-^, which 
provides the ratio [HCO+]/[CO] ~ 2 x 10"''. In the 
high velocity component, V > 570 km s""'^, the ratio is 
five times larger than in the low velocity components, 
[HCO+]/[CO] ~ 1 X 10-3. 

The [HCO+]/[CO] ratio is even larger within 
570 kms-i < F < 600 km s^^. While [HCO+]/[CO] 
~ 5 X 10"** for components at 1^ > 600 km s^^, in the 
velocity range between 570 km s^^ < ^ < 600 km s^^ 
it is a factor 10 larger, [HCO+]/[CO] ~ 5 x lO^^. This 
difference arises in the velocity range corresponding to 
the broad absorption component. 

5.2. Atomic versus molecular gas content 

In Figure 8 we present a comparison between the H I 
and C0(2-l) absorption features. It is important to note 
that we are tracing both the H I component d and molec- 
ular gas in the line of sight within a region located closeby 
(< 0.3 pc). 

A wide range of atomic to molecular abundance ra- 
tios are seen. While the main H I absorption features 
are located between V = 544 and 600 km s'^, in CO 
and other molecules the main absorption features occur 
close to the systemic velocity, and then at larger veloc- 
ities, V > 605 km s~^. The H I absorption features are 
remarkably more prominent at 560 < V^ < 610 km s^^ 
than their molecular counterparts. The two components 
5 and 6 of the C0(2-l) and HCO+(1-0) spectra at V 
> 605 km s^^ are located in the same velocity range as 
the H I broad absorption feature which is seen up to V 
= 640 km s-i. 

Next we compare the molecular to atomic abundance 
ratio assuming Tex = 10 K for C0(2-l) and Ts= 100 K 
for H I. First, in the low velocity component the gas is 
mostly in its molecular phase. The C0(2-l) line 1 seems 
to be related to the H I component IV, and it yields a 
ratio iV(H2)/^(H I) = (7.25 x lO^o cm'^) / (0.9 x 10^° 
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The A^(H2)/A^(H I) ratio corresponding to 



C0(2-l) line 2 is simfiar to that of the C0(2-l) line 
1. The H I absorption line SV and C0(2-l) line 3, the 
ratio of H2 to H I column densities is iV(H2)/iV(H I) > 



(1800 X 10 
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cm 



'■) I (12 X 10^ 



= 1.5. 



In the high velocity component located between 
570 km s-^< V < 600 km s^^, the gas is predominantly 
either atomic or molecular gas traced by HCO+(1-0), 
but not by C0(2-l). H I lines I/II and HI contain 
a total of N(R I) ~ (13.4 ± 2) x 10^° cm'^. If spa- 
tially linked to molecular clouds with the same kine- 
matics, they would only correspond to C0(2-l) line 4, 
which contains a smaller amount of molecular gas, A^(H2) 
~ (26 ± 4) X 10^^ cm-2. The ratio is iV(H2)/iV(H I) ~ 
0.02, 2-3 orders of magnitude smaller than in the low 
velocity components. However, we have to add the con- 
tribution from the molecular gas column density traced 
by other coexisting species such as HCO+(1-0) in this ve- 
locity range. Wiklind & Combes (1997) estimated it to 
be Af(HCO+)/A^(H I) ~ 1.5 x 10"^, or iV(H2)/A^(H I) 
~ 0.2 - 2 assuming a typical [HC0+]/[H2] = lO""^ - 
10^*. Even with the contribution of HCO+ in the ve- 
locity range between 570 km s-^< V < 600 km s~^, we 
obtain a slightly smaller [H2]/[H I] ratio than in the low 
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velocity component. 

Finally, from V > 600 km s^^, C0(2-l) components 
5 and 6 are within the H I broad line. There iV(H I) = 
7 X 10^" cm~^ (broad line), to be compared to the A^(Il2) 
~ 2.6 X 10^0 cm-^ which yields a ratio 7V(H2)/iV(H I) 
~ 0.4. 

6. DISCUSSION 

6.1. Warped Disk or Non-Circular / Infalling Motions? 

Although there is a consensus that narrow absorption 
lines with velocities close to the systemic velocity (Low 
Velocity Component -LVC-, in the velocity range 540 
km s~^< V < 560 km s~^, following the nomenclature 
by Wiklind & Combes 1997) represent gas far from the 
nucleus (~ 1 kpc), still there is no agreement regarding 
the more redshifted components (High Velocity Compo- 
nent -HVC-, 560 km s-^< V < 640 km s"!). 

A schematic view of proposed models to describe the 
almost edge-on disk of CenA is presented in Fig. 9. This 
figure includes edge-on and face-on views of the proposed 
models within the inner 2 kpc, as well as a qualitative 
interpretation of the kinematic location of the absorp- 
tion lines. A warped and thin disk model (e.g. Quillen 
et al. 2006; Neumayer et al. 2007; Quillen et al. 2009) 
reproduces well the observed emission lines in terms of 
distribution and kinematics. However, this model is not 
able to reproduce by itself the existence of the redshifted 
components (HVC) because it assumes pure circular or- 
bits (Eckart et al. 1999) (see Figure 9 a). It is necessary 
to augment this model with either diffuse gas (Eckart 
et al. 1999) (see Figure 9 b) or non-circular/infalling mo- 
tions (e.g., as suggested in van der Hulst et al. 1983 and 
Paper I; see Figure 9 c). Since the implications of the 
latter two models diverge, we discuss next what is the 
most plausible one. 

It is remarkable that we detect the broad line of the 
HVC against the brightest continuum source (position d) 
but it is found to be much weaker toward position c just 
~ 0.4 pc next to it and further from the nucleus. This 
considerably restricts the upper limit of a few hundred 
parsecs set by VLA data (Sarma et al. 2002). The non- 
circular/infalling motions within a warped disk model 
(Figure 9 c) can provide the observed differences between 
positions c and d, since one would expect higher veloc- 
ities closer to the nucleus. On the other hand, this is 
not expected in the warped disk with diffuse gas model 
(Figure 9 b) since this gas should be distributed at high 
latitudes up to about a few hundreds of parsecs. 

In addition, the physical locations of the different ab- 
sorbing regions in the line of sight derived from the 
warped disk with diffuse gas model (Figure 9 b) are com- 
plex. The model predicts that LVC lines 1 and 3 (Fig- 
ure 7) arise at r = 1.7 - 1.9 kpc and are separated by up 
to Az ~ 160 pc from the disk, line 2 should be located 
r = 200 pc and Az ~ 0, and the HVC lines 4i 5 and 6 
at r = 200 - 600 pc and Az ~ 300 pc (see Figure 9 6). 
In this model, within r < 600 pc, the nearer the corre- 
sponding absorbing region is to the nucleus the closer the 
corresponding line velocity is to the systemic velocity. 

Finally, no component within the inner r = 200 pc is 
taken into account in Eckart et al.'s model, although a 
large concentration of molecular gas is found there (Pa- 
per I) . A contribution to the molecular gas profile of this 
component should be seen as an extension of the narrow 



line 2 (arising from r — 200 pc at Az = pc in Eckart 
et al.'s model). The limit imposed by Eckart et al. (1999) 
was set at r = 200 pc so that the line 2 reproduces the 
observed one. However, without that constraint, line 2 
is difficult to be reproduced in that model. 

These arguments suggest that the nature of the HVCs 
is not just diffuse gas within a warped disk far from the 
nucleus, and that non-circular motions and/or infalling 
gas must be important mechanisms to explain the nature 
of this component. Within this scenario, molecular lines 
4, 5 and 6 and atomic lines I/II and HI would arise at 
radii r < 1 kpc, while lines 1, 2 and 3 and SV and IV 
would be at r > 1 kpc (Figure 9 c). 

6.2. An X-ray Dominated Region (XDR)? 

Since a powerful AGN resides in this galaxy, and a 
significant amount of gas is located within the inner r 
~ 220 pc (Paper I), it is reasonable to think that X-ray 
radiation plays a major role modifying the chemical com- 
plexity (i.e., dissociating molecular hydrogen into atomic 
gas, ionizing and heating, etc.). 

In fact, the AGN of Centaurus A is characterized by 
a X-ray luminosity of Lx ~ 4.8 x 10^^ erg s""'^. The 
incident flux at a distance of 1 pc is thus quite large, -Fx 
^ 300 erg s~^ cm~^. On the other hand, it is not likely 
that it has such a strong effect at distances ~100 pc, 
with fluxes of the order of Fx ^ 0.03 erg s~^ cm~^. The 
gas column density attenuating these X-rays is estimated 
to be ^ 10^'^ cm^^, probably by a cloud that entirely 
surrounds the nucleus or by a torus-like structure (Evans 
et al. 2004, and references therein). 

The LVC is likely arising from cold molecular gas, with 
kinetic temperatures T ^ 10 K and average hydrogen 
densities of the order of n ^ 10'' cm^'^, which seems 
to be appropriate for the disk as seen in emission (e.g. 
Eckart et al. 1990b) as well as in absorption (e.g. Wik- 
hnd & Combes 1997; MuUer & Dinh-V-Trung 2009). Less 
information is present in the literature for the physical 
conditions of the HVC, since this weak component is dif- 
ficult to detect. 

Observationally, Wiklind & Combes (1997) show that 
the average of [HCO+]/[HCN] - 1.3 and 1.7, while the 
average for [HNC]/[HCN] is 0.25 and 0.5, respective to 
the LVC and the HVC. The numerical ratios of Meijerink 
& Spaans (2005) for a XDR are [HCO+]/[HCN] - 1 - 10 
and [HNC]/[HCN] - 0.01 - 1, with Fx=160 erg s^^ cm^^ 
and uh — 10^'^ cm~'^ (or both Fx and uh two orders 
of magnitude lower). This numerical ratio agrees well 
with the observations. However, under certain circum- 
stances, a PDR could also reproduce the given molecular 
gas ratios. 

X-rays do not lead to strong dissociation of CO and 
it can be present in an XDR at elevated temperatures 
(Spaans 2008). Warm CO gas produces emission origi- 
nating from high rotational transitions. Bright CO emis- 
sion lines have been seen for the first three J levels, and 
larger T^x (20 - 30 K) have been inferred for molecular 
clouds close to the nucleus (e.g. Israel et al. 1991). On 
the other hand, absorption lines of the lower J transitions 
will preferentially sample the (excitationally) coldest gas. 
In the nuclear regions, where higher densities and kinetic 
temperatures are expected, only high- J transitions of CO 
would have such broad lines as seen in HCO+(1-0), as 
it has been suggested with previous CO (3-2) single-dish 
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data (Wiklind & Combes 1997). This is likely the reason 
why we did not detect the broad line in C0(2-l), but 
may well be detected in higher transitions of CO. This 
could be a good test to distinguish the physical condition 
of the molecular gas closer to the nucleus. 

The molecular lines at V^ > 605 km s~^, including ab- 
sorption features 5 and 6, are probably molecular clumps 
close to or within the H I region from which the H I broad 
absorption component arises. We consider the possibil- 
ity that these features are located close to the nucleus 
(as suggested in Figure 9 c). For an X-ray illuminated 
gas like these clouds, the physical conditions of the gas 
that surrounds the AGN can be studied using the effec- 
tive ionization parameter, ^off, which is proportional to 
the ratio of incident (and attenuated) X-ray photon flux 
to gas density of the cloud (Maloney et al. 1996). The 
effective ionization parameter of a molecular cloud can 
be calculated as: 

Coff - 1.1 X IQ-^LM/iN^i'ngr^) 

where the L44 = ix/lO^'* is the X-ray luminosity of hard 
photons (energies > 2 keV) in ergs s~^, ng the gas den- 
sity of the cloud in 10^ cm^'^, r the distance to the X-ray 
emitting source in pc, and N22 the attenuating gas col- 
umn density in 10^^ cm^^ (Maloney et al. 1994, 1996). 
The molecular cloud regions are exposed to an X-ray lu- 
minosity from Cen A's AGN of Lx = 4.8 x 10"*^ ergs s~^ 
in the hard energy range 2 -10 keV (Evans et al. 2004). 

If we assume for the molecular gas arising in the two 
high velocity absorption features 5 and 6 a Tex — 10 K, 
a size of the 1.3 mm continuum of 0.01 pc (Kellcrmann 
et al. 1997), and a size of the cloud of R ~ 10 pc then 
the molecular clouds would be characterized by densities 
of about n ^ 10^ cm~^. 

The gas column densities that arc attenuating the X- 
rays are estimated to be ~ 10^"^ cm^^ (Evans et al. 2004). 
By assuming a distance to the nucleus of r = 10 pc, then 
it yields a parameter ^eff = 7 xlO"**. For r = 1 pc then 
^eff = 0.1. A ^off > 10""^ means that the gas will be pre- 
dominantly in its atomic phase and it indicates that the 
molecular clouds must be at r > 10 pc in order to sur- 
vive under the above assumptions. This is in agreement 
with the smaller [H2]/[H I] ratio found in the HVC with 
respect to the LVC (§ 5.2). 

6.3. Physical properties of the circumnuclear HI 

As mentioned in § 3.2 the velocity widths of our three 
main prominent H I features (i.e. I, I/II and III) are 
between 5-10 km s"'^. If one neglects line broadening 
due to rotation and turbulence, and attributes the ob- 
served FWHM only to thermal motions in the gas, the 
line widths shown in Table 2 would indicate tempera- 
tures between 800 and 5,700 K. Turbulence can account 
for a linewidth of about ~ 5 km s~^, at least within our 
own Galaxy (Burton 1988). Taking this into account, 
the actual gas temperature is probably between 20 and 
2500 K. 

The nature of the gas from which the broad absorp- 
tion lines arise are much less understood. The underlying 
broader H I absorption line, with a width of '-^ 55 km s"^, 
cannot be attributed solely to thermal motions since it 
would give too high temperatures (~ 60,000 K). As in- 
dicated by Sarma et al. (2002) this redshifted broad H I 



component is probably a blend arising from several com- 
plexes rather than a single one. Since it only shows 
redshifted velocities, it is likely that this component is 
largely affected by systematic kinematic contributions, 
such as non-circular and/or infalling motions, and close 
to the nucleus, mostly within a distance < 0.7 pc (0.4 pc 
from position c to d, and 0.3 pc from position d to the 
actual nucleus; as mentioned in § 3.1). If the geometry 
of the absorbing atomic gas is similar to the molecu- 
lar gas disk/torus in Paper I (inclination i = 70°), then 
the disk-like feature associated with the broad H I line 
with column densities of about iV(H I) ~ 7 x 10^° cm^^ 
would have a (non-projected) radius r < 2 pc. However, 
we cannot rule out other geometries. Also, note that 
the kinematic coincidence between the broad absorption 
lines of H I and other molecular lines such as HCO+(1-0) 
may suggest that the corresponding absorbing molecular 
regions are physically connected. 

The atomic material arising from the broad compo- 
nent is likely contributing to fuel the powerful AGN of 
Centaurus A. Assuming that this H I region is infalling 
with a velocity of ~ 37 km s~^(mean velocity minus the 
systemic velocity), a column density of A^(II I) ^ 7 x 
10^° cm~^, and a size of 2 pc, we estimate an accretion 
rate of 2 x 10~^ M0 yr~^, enough to power the required 
accretion rate of ~ 1.5 x 10~^ M© yr~^ obtained from its 
isotropic radio luminosity (e.g. van Gorkom et al. 1989). 

7. SUMMARY AND CONCLUSIONS 

We present H I and C0(2-l) absorption features to- 
ward the nuclear regions of Centaurus A (NGC 5128), 
using the VLBA and the SMA respectively. 

Our H I observations allow us to discern with sub-pc 
resolution the absorption profiles toward the nuclear jet 
in the inner 5.4 pc. The H I absorption lines are com- 
posed of a system of 4 main narrow lines as well as an 
underlying broad line (FWHM ~ 55 km s~^), in agree- 
ment with previous works using lower spatial resolution 
observations. Interestingly, the underlying broad absorp- 
tion line in our data is seen to be more prominent toward 
the central 21 cm continuum component (base of the nu- 
clear jet, although not the nucleus itself) than toward 
the continuum emission at > 0.4 pc from it and further 
from the nucleus. 

We interpret this as a broad velocity component that 
arises from an absorbing region that is physically close 
to the AGN itself, and not just as diffuse gas seen in 
projection as suggested by Eckart et al. (1999) since the 
height of such diffuse component is expected to be of the 
order of a few hundreds of parsec, and thus we would not 
expect any substantial difference between both positions 
at a 0.4 pc scale. In addition, the large width of the 
broad absorption line suggests that the absorbing region 
is composed of several complexes with systematic non- 
circular or infalling motions. 

As for the molecular gas component, our C0(2-l) in- 
tcrferometric data provides for the first time CO absorp- 
tion features against the unresolved 1.3 mm continuum 
emission with minimal contamination by line emission. 
The narrow lines (likely arising from clouds at large radii 
from the nucleus) observed in the C0(2-l) and ^'^C0(2- 
1) profile were identified with previously observed molec- 
ular lines such as HCO+(1-0) (e.g. Wiklind & Combes 
1997; van Langevelde et al. 1995). Overall, our calculated 
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isotopic ratios [^'^C]/[^^C] ~ 5 - 17 and the [^^0]/[^^0] 
~ 70 - 220 seem to indicate that the molecular clouds 
producing the observed narrow absorption lines close to 
the systemic velocity have been in a relatively long-term 
starburst environment as in nearby starburst galaxies. 
We did not detect any counterpart of the broad absorp- 
tion line in the CO (2-1) spectrum as that found in H I 
or other molecular lines, although the noise level in our 
C0(2-l) profile would have been sufficient to detect this 
feature if it were as strong as that of HCO+(1-0). 

Since both our H I (position d) and molecular absorp- 
tion profiles arise from molecular clouds in front of a sim- 
ilar continuum emission region within ~ 0.3 pc, we can 
compare the atomic and molecular gas properties along 
the line of sight. The abundances of H I to H2 are no- 
ticeably different at different velocities. While the low 
velocity components, V < 560 km s~^ are dominated by 
molecular gas, the atomic phase gas seems to dominate 
at 560 < V < 600 km s~^ when the H I line is compared 
with the molecular gas content as traced by the C0(2- 
1) line. Although the molecular to atomic gas ratio is 
not drastically different when we consider the molecu- 
lar gas traced by IICO+(1-0), still there is a trend for 
a larger atomic content for components at higher veloci- 
ties. These distinct signatures suggest that the physical 
properties, location and chemistry of the gas producing 
the broad line is different to those of the narrow lines. 

Since the gas corresponding to the broad line is likely 
not far from the AGN, we argue that its properties might 
be modified by X-rays, which penetrate much farther 
than UV. The observed abundances of species such as 
IICO+ and HCN were estimated to be ~ 10"^ , and 
abundance ratios [HCO+]/[HCN] > 1 and [HNC]/[HCN] 
~ 0.5, which seems to agree with the chemistry of a X-ray 



dominated region (XDR). In the nuclear regions, higher 
densities and kinetic temperatures are expected. Only 
high-J transitions would have such broad lines as seen 
in IICO+(1-0). Since X-rays do not strongly dissociate 
CO molecules and they can coexist at high temperatures, 
higher transitions of CO other than C0(2-l) could trace 
the broad line component, as has been suggested with 
previous CO (3-2) single-dish data. 

Because of its kinematics, the two narrow red-shifted 
lines at V^ > 600 km s^^ found in our C0(2-l) spectrum, 
as well as in previous HCO'*'(1-0) spectrum might be 
close to the AGN and falling toward it. These molecular 
clouds are likely at a distance > 10 pc from the AGN so 
that they can survive the X-ray emission. On the other 
hand, at distances < 10 pc most of the gas will be in 
its atomic phase, material that is likely contributing to 
power the AGN of Centaurus A. 



We thank the SMA and NRAO staff members who 
made the observations reported here possible. We also 
thank A. Sarma for providing the 21 cm continuum VLA 
data, and S. Martin and I. Jimenez-Serra for interest- 
ing discussions. This research has made use of NASA's 
Astrophysics Data System Bibliographic Services, and 
has also made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space 
Administration. DE was supported by a Marie Curie 
International Fellowship within the 6**^ European Com- 
munity Framework Programme (MOIF-CT-2006-40298). 

Facilities: VLBA,SMA 



REFERENCES 



Burton, W. B. 1988, The structure of our Galaxy derived from 

observations of neutral hydrogen (Galactic and Extragalactic 

Radio Astronomy), 295—358 
Combes, F. 2008, Ap&SS, 313, 321 
Eckart, A. et al. 1990a, ApJ, 365, 522 
— . 1990b, ApJ, 363, 451 

Eckart, A., Wild, W., & Ageorgcs, N. 1999, ApJ, 516, 769 
Espada, D. ct al. 2009, ApJ, 695, 116 
Evans, D. A. ct al. 2004, ApJ, 612, 786 
Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 

590 
Gardner, F. F., & Whiteoak, J. B. 1976a, Proceedings of the 

Astronomical Society of Australia, 3, 63 
— . 1976b, MNRAS, 175, 9P 
Gupta, N. et al. 2007, ApJ, 654, Llll 
Harris, G. L. H., Rejkuba, M., & Harris, W. E. 2009, ArXiv c-prints 

0911.3180 
Harrison, A., Henkcl, C, & Russell, A. 1999, MNRAS, 303, 157 
Henkel, C. et al. 1993, A&A, 268, L17 
— . 2009, A&A, 500, 725 

Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 
Horiuchi, S., Meier, D. L., Preston, R. A., & Tingay, S. J. 2006, 

PASJ, 58, 211 
Israel, F. P. 1998, A&A Rev., 8, 237 
Israel, F. P., Raban, D., Booth, R. S., & Rantakyro, F. T. 2008, 

A&A, 483, 741 
Israel, F. P., van Dishoeck, E. F., Baas, F., de Graauw, T., & 

Phillips, T. G. 1991, A&A, 245, L13 
Israel, F. P. et al. 1990, A&A, 227, 342 
Jones, D. L. et al. 1996, ApJ, 466, L63 
Kellermann, K. I., Zensus, J. A., & Cohen, M. H. 1997, ApJ, 475, 

L93 
Maloney, P. R., Bcgclman, M. C, & Rees, M. J. 1994, ApJ, 432, 

606 



Maloney, P. R., HoUcnbach, D. J., & Tielens, A. G. G. M. 1996, 

ApJ, 466, 561 
Meijerink, R., & Spaans, M. 2005, A&A, 436, 397 
Meisenheimer, K. et al. 2007, A&A, 471, 453 
MuUer, S., & Dinh-V-Trung. 2009, ApJ, 696, 176 
MuUer, S., GueUn, M., Dumke, M., Lucas, R., & Combes, F. 2006, 

A&A, 458, 417 
Neumayer, N. et al. 2007, ApJ, 671, 1329 

Prantzos, N., Aubcrt, O., & Audouze, J. 1996, A&A, 309, 760 
Quillcn, A. C. et al. 2008, MNRAS, 384, 1469 
— . 2006, ApJ, 645, 1092 
Quillcn, A. C, Neumayer, N., Oosterloo, T., & Espada, D. 2009, 

ArXiv e-prints 0912.0632 
Roberts, M. S. 1970, ApJ, 161, L9 

Sarma, A. P., Troland, T. H., & Rupen, M. P. 2002, ApJ, 564, 696 
Schiminovich, D., van Gorkom, J. H., van dcr Hulst, J. M., & 

Kasow, S. 1994, ApJ, 423, LlOl 
Seaquist, E. R., & Bell, M. B. 1990, ApJ, 364, 94 
Shepherd, M. C., Pearson, T. J., & Taylor, G. B. 1995, in Bull. 

Amer. Astron. Soc, Vol. 27, , 903 
Solomon, P. M., Sanders, D. B., & Scoville, N. Z. 1979, ApJ, 232, 

L89 
Spaans, M. 2008, EAS Publications Series, Vol. 31, cd. C. Kramer, 

S. Aalto, & R. Simon, 47-520 
Tingay, S. J. et al. 1998, AJ, 115, 960 
Tingay, S. J., & Murphy, D. W. 2001, ApJ, 546, 210 
van der Hulst, J. M., Golisch, W. F., & Haschick, A. D. 1983, ApJ, 

264, L37 
van Gorkom, J. H. et al. 1989, AJ, 97, 708 
van Gorkom, J. H., van der Hulst, J. M., Haschick, A. D., & Tubbs, 

A. D. 1990, AJ, 99, 1781 
van Langevelde, H. J., Pihlstrom, Y., & Beasley, A. 2005, Ap&SS, 

295, 249 



10 



D. Espada et al. 



^ 0.5 



<u 

N 

•rH 

o 



0.5 



tl{ > ^w\f -v^*' 



^^■VTM 



1 — \ — M — I— 



yi 



' ^■^VMIK^V^"'"^ 



a) 

All visibilities 



»vtiftV^ff^iJLfyy^' 



- 12 



co(2-i; 



^'*'ipf^liifrH^^'^^""tf'* 



b) 

Visibilities > 50m 



_i \ i_ 



_i 1 \ i_ 



_L 



400 500 600 700 

LSR Velocity (km s"^) 



Fig. 1. — a) CO(2-l) normalized intensity spectra obtained from all the visibilities using a fit in the UV plane assuming a point-like 
feature for the continuum emission. It clearly fails to eliminate most of the emission arising from the circumnuclear molecular gas (Espada 
et al. 2009). The arrow indicates the systemic velocity of the galaxy (V = 541.6 km s~^). b) Same as a) but only using visibilities with 
projected baselines larger than 50 m. The contamination by emission is efficiently diminished. 
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Fig. 3. — H I absorption profile of Cen A, integrated over the entire 21 cm continuum region probed by our VLBA experiment (Fig. 2). 
The profile is mainly composed of five lines: IV, SV, I/II and III, plus an underlying broader line with a width of 53 ± 10 km s~^, whose 
wing is clearly seen at V > 600 km s~^. The velocity resolution is 3.3 km s~^ and the rms noise of 3 mjy beam~^ per channel. The solid 
line is the sum of the Gaussian fits to the individual lines (dashed lines, as specified in Table 2). The plot at the bottom shows the residual 
from this fit, with a rms a = 0.01 Jy. 
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Fig. 4. — H I absorption lines toward single pixels at positions a, b, c and d as outlined in Figure 2 (see bottom right of each panel), as 
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TABLE 1 
Main parameters of the VLBA H I and SMA C0(2-1) observations 



H 1 VLBA observations 



Date 

R.A. of phase center (J2000) 
Decl. of phase center (J2000) 
Time on source (hr) 
FWHM of synthesized beam 



1994 November 26 

13'^25™27.=609 

-43°01'07y97 

2 

36.3 X 7.2 mas (0.65 x 0.13 pc), PA = 14?8 



LSR velocity at band center (km s ^) 541.6 

Total bandwidth 2 MHz (~ 430 km s^^) 

Spectral resolution (km s^-"^) 3.3 

rms noise (3.3 km s^^) 3 mjy beam~^ 

rms (continuum map) 2 mJy beain"^ 

CO(2-l) SMA observations 



Date 

R.A. of phase center (J2000) 

Decl. of phase center (J2000) 

On source time (hr) 

FWHM of synthesized beam 

LSR velocity at band center (km s^-*^) 

Total bandwidth 

Spectral resolution (km s^^) 
rms noise (1 km s^^) 



2006 April 5 

13'^25™27.=6 

-43°01'09" 

3 

2('4 X 6^'0 (43 X 108 pc), PA = 0?82 

550 

2 GHz (~ 2870 km s^^) in each sideband 

(separated by 10 GHz) 

1 

0.1 Jy bcam^-'^ (0.3 Jy bcam^-*^ for baselines > 50 m) 



TABLE 2 
Parameters of H I absorption lines integrated over the entire 21 cm continuum 



Peak S ('') 
(mJy beam^^) 



■pWHMf^ 
(km s^-"-) 



Absorption V 

feature (km s"-*^) 



Peak T 



JrdV 
(km s -*-) 



Af(HI)/T, 
(cm-2 K-i) 



IV 


541 ± 1 


-57 ± 19 


sv 


552 ± 1 


-323 ± 20 


I, II 


575 ± 1 


-203 ± 17 


HI 


595 ± 1 


-122 ±20 


Broad 


578 ± 4 


-34 ±21 



5 ± 2 0.13 ± 0.03 0.7 ± 0.4 

8.5 ± 1.2 0.74 ± 0.03 7.7 ± 1.0 

12.6 ± 1.6 0.48 ± 0.02 7.6 ± 0.9 

5.2 ± 1.4 0.25 ± 0.03 1.8 ± 0.4 

53 ± 10 0.07 ± 0.03 4.0 ± 1.9 



(1.2 ±0.7) X 10^* 
(14 ±2) X 10^** 
(14 ±2) X 10^** 

(3.2 ±0.7) X lO^* 
(7 ±3) X 10^** 



^Flux density peak and FWHM are calculated from a Gaussian fit. 
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TABLE 3 

Parameters of H I absorption lines toward individual positions ('') 



Absorption 
feature 



V 
(km s^^) 



Peak S 
(mjy beaiii^^) 



FWHM 
(km s~^) 



Peak T 



(km s ^) 



N{n l)/T, 

(cm-2 K-i) 



a) 
IV 

I/II 

III 

Broad 



(541) 
552 ± 1 
(575) 
(595) 
(578) 



< 2 
~2±3 

< 2 

< 2 

< 2 



(5) 

4 

(14) 

(6) 

(50) 



< 0.6 

< 0.6 

< 0.6 

< 0.6 

< 0.6 



< 2 

< 2 

< 3 

< 2 

< 6 



< 4 X 10^* 

< 4 X 10^** 

< 11 X 10i« 

< 5 X 10^*^ 

< 43 X 10^* 



b) 














IV 


(541) 


< 17 


(5) 


< 0.3 


< 1 


< 2 X 10^* 


sv 


552 ± 1 


9±2 


5 ± 3 


0.81 ± 0.07 


4 ± 2 


(7 ± 3) X 10^** 


I/II 


(575) 


< 17 


(14) 


< 0.3 


< 2 


< 7 X 10^** 


III 


(595) 


< 17 


(6) 


< 0.3 


< 1 


< 4 X 10^* 


Broad 


(578) 


< 17 


(50) 


< 0.3 


< 4 


< 30 X 10^** 



IV 


541 ± 1 


129 ±5 


5 ± 2 


0.12 ± 0.02 


0.7 ± 0.3 


(1.3 ± 0.5) 


X 


10i» 


SV 


551 ± 1 


59 ±5 


6.3 ± 1.0 


0.90 ± 0.02 


5.5 ± 1.0 


(10 ± 2) 


X 


10i« 


I/II 


576 ± 1 


108 ±5 


13 ± 2 


0.29 ± 0.02 


4.4 ± 0.7 


(8.0 ± 1.3) 


X 


10i» 


III 


596 ± 1 


120 ±5 


5 ± 2 


0.18 ± 0.02 


1.2 ± 0.4 


(2.2 ± 0.7) 


X 


10" 


Broad 


578 ±4 


142 ±5 


53 ± 10 


0.01 ± 0.02 


0.5 ± 1.1 


(1±2) 


X 


10i» 


d) 
IV 


540 ± 1 


546 ± 10 


3 ± 2 


0.10 ± 0.03 


0.5 ± 0.2 


(0.9 ± 0.4) 


X 


1018 


SV 


551 ± 1 


286 ± 10 


7.7 ± 1.3 


0.75 ± 0.03 


6.7 ± 1.1 


(12 ± 2) 


X 


1018 


I/II 


574 ± 1 


400 ± 10 


12.8 ± 1.8 


0.41 ± 0.03 


5.9 ± 0.9 


(10.7 ± 1.7) 


X 


1018 


III 


595 ± 1 


506 ± 10 


5.6 ± 1.6 


0.18 ± 0.03 


1.5 ± 0.4 


(2.7 ± 0.7) 


X 


1018 


Broad 


577 ±4 


564 ± 10 


55 ± 10 


0.06 ± 0.03 


3.6 ± 1.3 


(7 ±2) 


X 


1018 



'^ Velocity {V), flux density peak (Peak S), FWHM (fufl width half maximum), peak optical depth (Peak t) and integrated optical depth 
(J T dV) were derived from Gaussian fits. Velocities in parenthesis are assumed values in the fitting. Uncertainties in the flux density peak 
and the peak optical depth correspond to the residuals after Gaussian fitting (see Fig. 4). Uncertainties in the integrated optical depth are 
calculated as the uncertainty of the area of the Gaussian profile. Upper limits in the integrated optical depth are 3a-, and are calculated as 
3 X cr X (Ay (5v)i'^, where AV is the velocity of the integrated profile and &v the channel spacing. Line SV in position a is detected, but 
the optical depth cannot be calculated since the absorption depth is larger than the continuum level due to the noise. Thus, we consider 
it as an upper limit. 
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TABLE 4 
Parameters of the C0(2-1), ^^C0(2-1) and C-'-*0(2-1) absorption lines 



Isotope 


Line 


V 
(km s-l) 


S 
(Jy beam~^) 


Peak T 


JrdV 
(kms-i) 


Af(CO)(") 
(1014 cm-2) 


Af(H2) (') 
(lO^* cm-2) 


12C0{2-1) 


1 


539 


0.72 ± 0.25 


2.1 ± 0.3 


6.00 ± 0.10 


616 ± 10 


725 ± 12 




2 


544 


1.20 ± 0.25 


1.6 ± 0.2 


2.33 ± 0.07 


239 ± 7 


281 ± 8 




3 


550 


0.16 ± 0.25 


> 3.6 


> 14.88 


> 1530 


> 1800 




4 


575 


5.33 ± 0.25 


0.10 ± 0.05 


0.22 ± 0.03 


23 ± 3 


26 ± 4 




5 


606 


2.66 ± 0.25 


0.80 ± 0.10 


1.68 ± 0.07 


172 ± 7 


203 ± 8 




6 


613 


4.61 ± 0.25 


0.25 ± 0.06 


0.48 ± 0.07 


49 ± 7 


58 ± 8 


i^CO{2-l) 


1 


541 


5.30 ± 0.10 


0.11 ± 0.02 


0.35 ± 0.03 


36 ± 3 


3600 ± 300 




2 


545 


4.90 ± 0.10 


0.19 ± 0.02 


0.46 ± 0.03 


47 ± 3 


4700 ± 300 




3 


552 


3.33 ± 0.10 


0.58 ± 0.03 


1.42 ± 0.04 


145 ± 4 


14500 ± 400 




4 


- 


- 


- 


<0.08 


<8 


<800 




5 


607 


5.80 ± 0.10 


0.05 ± 0.01 


0.10 ± 0.02 


10 ± 2 


1000 ± 200 




6 


- 


- 


- 


<0.09 


<9 


<900 


Ci«0(2-1) 










<0.11 


< 11 


< 6600 



Tsotopic CO column densities. We assume a filling factor of unity, local thermodynamic equilibrium (LTE) and an excitation temperature 
Tex = 10 K for the absorbing molecular gas (§ 4.2). Note that the CO(2-l) line 3 is likely saturated, so the optical depth and the column 
densities are lower limits. ^^CO lines 4 and 6 and C^^O lines are not detected, thus we use upper limits taken at a 3 ct level. 

''H2 column density (Af(H2)) derived from the CO(2-l) absorption lines (and isotopes). To derive the A'^(H2) from the CO(2-l) lines, we 
use an abundance ratio of [CO]/[H2] = 8.5 X IQ-^ (Frerking et al. 1982), [^^CO]/[U2] = 1 X 10"'' (Solomon et al. 1979) and [Ci*0]/[H2] 
= 1.7 X 10-7 (Frerking et al. 1982) ([CO]/[i3cO] = 85 and [CO]/[Ci«0] = 500]). 



